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Abstract: We study the question of diagonalizability of the Hamiltonian for the Faddeev- 
Reshetikhin (FR) model in the two particle sector. Although the two particle S-matrix ele- 
ment for the FR model, which may be relevant for the quantization of strings on AdS^ x S^, 
has been calculated recently using field theoretic methods, we find that the Hamiltonian 
' for the system in this sector is not diagonalizable. We trace the difficulty to the fact that 

■ the interaction term in the Hamiltonian violating Lorentz invariance leads to discontinuity 

■ conditions (matching conditions) that cannot be satisfied. We determine the most gen- 
. eral quartic interaction Hamiltonian that can be diagonalized. This includes the bosonic 

Thirring model as well as the bosonic chiral Gross-Neveu model which we find share the 
^ I same S-matrix. We explain this by showing, through a Fierz transformation, that these 

H ' two models are in fact equivalent. In addition, we find a general quartic interaction Hamil- 

tonian, violating Lorentz invariance, that can be diagonalized with the same two particle 
S-matrix element as calculated by Klose and Zarembo for the FR model. This family of 
generalized interaction Hamiltonians is not Hermitian, but is PT symmetric. We show 
that the wave functions for this system are also PT symmetric. Thus, the theory is in 
a PT unbroken phase which guarantees the reality of the energy spectrum as well as the 
unitarity of the S-matrix. 



Keywords: Sigma Models, AdS-CFT and dS-CFT Correspondence, Integrable Field 



Theories 



*Also at Saha Institute for Nuclear Physics, 1/AF Bidhannagar, Calcutta 700064, India 



Contents 

|l]. Introduction ^ 

|2|. The Faddeev-Reshetikhin model ^ 

H. String Sigma Model. g 

^ Diagonalization of the Hamiltonian [l^ 

|5|. Solving the Boundary Condition 17 

|6|. Diagonalizability and PT Symmetry 21 

|7|. Conclusion 25 



1. Introduction 

The classical integrability of the superstring on AdS^ x [1"8] has led to a lot of in- 
teresting studies recently. Through AdS/CFT correspondence, a lot is already known on 
the gauge theory side (for reviews see [9-17] and references therein). However, the string 
theory presents several technical difficulties, as a result of which, even though it is believed 
that integrability should hold in the quantum theory, quantizing the string remains an 
open question until now. It is worth recalling here that the Green-Schwarz string can be 
described by a symmetric space sigma model. The flat currents in this model [1, 18, 19] 
which define the basic variables of the theory, like all sigma models, satisfy non-ultralocal 
Poisson bracket structures (that involve derivatives of delta functions) [20-24]. This is one 
of the main difficulties in carrying out the quantization of this model [22,25-31]. Quanti- 
zation of the model is, of course, absolutely crucial in understanding, say, the spectrum of 
the theory. And this still remains an open question. 

In the context of the principal chiral sigma model on SU{2), Faddeev and Reshetikhin 
[32] have suggested that the question of quantization may be carried out in the following 
manner. They propose that the Hamiltonian for the original sigma model as well as the non- 
ultralocal Poisson brackets may be replaced by another Hamiltonian and a new Poisson 
bracket structure (that is ultralocal) which lead to the same dynamical equations. One 
should quantize this new system and then recover the original sigma model afterwards in 
a certain limit. The idea is completely parallel to that contained in the method of Bethe 
ansatz in relativistic models (such as the massive Thirring model), where one studies the 
S-matrix elements by expanding the theory around the wrong vacuum and then tries to go 
back to the true vacuum of the theory where all the negative energy states are filled. The 
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simpler (FR) model proposed by Faddeev and Rcshetikhin is quite important from this 
point of view in connection with quantization of strings on AdS^ x [2,28,31,33-40]. 

Recently, Klose and Zarembo (KZ) [41] calculated the S-matrix for a number of 1 + 1 
dimensional integrable models (see also [42-46]) using standard field theoretic methods 
[47-50] in a simple manner. The simplicity of their method arises from the fact that the 
calculations are carried out in the wrong vacuum [51,52]. In this case, it is well known, 
for example, that in the two particle sector, the contribution to the S-matrix comes only 
from the bubble diagrams and if the system is integrable, all other scattering elements 
can be related to the two particle S-matrix [53-56]. One of the models studied by KZ is 
indeed the FR model, whose S-matrix element for the positive energy two particle states 
has a simple form that reflects the violation of Lorentz invariance present in the interaction 
Hamiltonian. This is interesting and, in fact, is relevant as a first step in understanding 
the quantization of the string itself. However, since the S-matrix element for the FR model 
is calculated in the wrong vacuum, it is necessary, as a next step, to go to the true vacuum 
to extract physical results [57]. This can be carried out by diagonalizing the (quartic) 
Hamiltonian of the theory in this sector. 

In this paper, we study the question of diagonalization of the two particle Hamiltonian 
for the FR model systematically. Surprisingly, we find that the quartic Hamiltonian for 
the FR model cannot be diagonalizcd. The problem arises because of the Lorentz violating 
term in the interaction Hamiltonian which leads to boundary (matching) conditions that 
cannot be satisfied to determine the wave function. This is rather puzzling given the 
nice S-matrix result of KZ. However, our results are consistent with the results of KZ in 
the following way. We find that, while the discontinuity cannot be matched across the 
boundary to determine the wave function for the system, the extra term violating the 
boundary condition is orthogonal to the positive energy two particle state. Therefore, it 
drops out when the inner product with positive energy states is taken and the discontinuity 
condition, in this case, yields the correct S-matrix obtained by KZ. This is completely 
consistent with the field theoretic calculation of KZ, which involves only a calculation of 
matrix elements. This, therefore, leads to the first important new feature that results 
from our analysis and which had not been observed earlier in other integrable models. 
Namely, while diagonalization of a system (Hamiltonian) leads to the S-matrix, having 
the S-matrix element (say, from a field theoretic calculation) does not automatically imply 
diagonalizability of the system. 

As a result of this lack of diagonalizability of the quartic Hamiltonian, we then searched 
for and determined the most general quartic Hamiltonian in this context for which the 
boundary conditions can be matched. The set of potentials which can thus be diagonalized 
include the bosonic Thirring model [58, 59] (which was known earlier to be integrable) as 
well as the bosonic chiral Gross-Neveu model (which to the best of our knowledge had 
not been studied earlier). Both these models respect Lorentz invariance and, in fact, we 
find that the two systems share the same S-matrix. Following this puzzling coincidence, we 
investigate and show, through a Fierz transformation, that the bosonic Thirring model and 
the bosonic chiral Gross-Neveu model are in fact equivalent. In addition, we determine a 
general quartic Hamiltonian, violating Lorentz invariance, which can be diagonalized and 
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leads to the same S-matrix as calculated by KZ for the FR model. We emphasize here 
that a field theoretic calculation with this generalized potential (interaction vertex) would 
yield the same S-matrix element as for the FR model. This is indeed quite interesting and 
another important result of our analysis, namely, different potentials (interaction vertices) 
can lead to the same S-matrix element in a field theoretic calculation. Furthermore, we 
find that in spite of the fact that the spectrum of this generalized family of Hamiltonians 
is real and the S-matrix is unitary, the Hamiltonian is not Hermitian. On closer analysis, 
we find that the new family of Hamiltonians is, in fact, PT symmetric [60-62] (for reviews 
see [63,64] and references therein). We show that the theory is in the PT unbroken phase 
which guarantees the reality of the spectrum as well as the unitarity of the S-matrix. This 
is the third important result of our analysis which identifies the relevance of PT symmetry 
with this integrable system. This, of course, still leaves us with the interesting question of 
why the FR model cannot be diagonalized and this is presently under study. 

Our paper is organized as follows. In section 2, we give a brief review of the Faddeev 
and Reshetikhin model and its relevance to the regularization of the ambiguities in the 
current algebra of the SU{2) principal chiral model. In section 3, we recapitulate the 
relation of the FR model to the string sigma model, in particular, within the context of 
the Rx subsector of AdS^ x background for simplicity. Here we also describe briefly 
the field theoretic calculation by KZ of the S-matrix for this model. In section 4, we 
demonstrate the non-diagonalizability of the corresponding Hamiltonian in the operator 
formalism and try to make connection with the field theoretic calculation. In section 5, we 
analyze the underlying quantum mechanical system to understand the difficulty in more 
detail. We show that the term in the interaction Hamiltonian, violating Lorentz invariance, 
leads to matching conditions that cannot be satisfied. Here we also make connection with 
the field theoretic results. In section 6, we present the general quartic Hamiltonian that 
can be diagonalized and determine the S-matrix associated with this system. Various 
special cases are studied here. In particular, we show that the bosonic Thirring interaction 
is equivalent to the bosonic chiral Gross-Neveu interaction. In addition, we determine a 
general Hamiltonian, violating Lorentz invariance, that can be diagonalized and has the 
same S-matrix as that calculated by KZ for the FR model. We show that this family 
of generalized Hamiltonians is not Hermitian, but is PT symmetric. By studying the 
transformation properties of the wave functions of the system, we show that the system 
is in the PT unbroken phase which guarantees the reality of the spectrum as well as the 
unitarity of the S-matrix. In section 7, we give a brief summary of our results. 

2. The Faddeev- Reshetikhin model 

The Lagrangian for the SU{N) principal chiral model is written in the form: 



where = -g-^d^g = J'^t" , [t«,t^] = o, 6, c = 1, 2, • • • A^^ - 1, c, g SU{N). Here 



7 is a constant and throughout the paper we use the Bjorken-Drell metric which, in 1 -|- 1 
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dimensions, has the diagonal form rj^'^ = (+,—). From the definition, we see that the 
current is a pure gauge and, therefore, satisfies a zero curvature condition in addition to 
the dynamical equation. Thus, it can be shown that the (current) variables of the theory 
satisfy the equations: 

5^ J'^''^ = 0, 

d,j:-d,r^ + r''jf,r, = o. (2.1) 

We can choose Jf to be the dynamical variable and carry out the Hamiltonian analysis 
for the system which leads to the following Poisson bracket structures for the dynamical 
variables of the theory: 

{JS{x),4iy)}=jr'^J^{x)6ix-y) 

{JSix), Jliy)} = Tr'Vf (x)5(x -y)- j5''%6ix - y) (2.2) 
W(x),Jt(y)} = 

Because of the presence of terms involving derivatives of delta functions, the Poisson bracket 
algebra becomes non-ultralocal. (We remark here parenthetically that the Poisson bracket 
for the nonlinear sigma model is not unique [65]. However, in any form, there always exist 
some derivatives of delta function, making them non-ultralocal.) This leads to ambiguities 
in the calculation of the basic algebra of transition matrices and renders inapplicable the 
standard procedure of quantization for such systems. 

The Faddeev-Reshetikhin (FR) model studies the quantization of the SU{2) principal 
chiral model for which a,b,c = 1,2,3, f"'^'^ = —e"-^'^ and the generators are related to the 
Pauli matrices, namely, = ia""/!. In this case, the basic variables of the theory can also 
be written as (3 dimensional) vectors in the internal space. Thus, defining 

5+ = i^(Jo + Ji), 



S. = i^(Jo - Ji 



(2.3) 



we note that in these variables, the equations in ( p.lj ) take the forms 

dtS± ^ d^S± ± 2jS+ xS-=0, (2.4) 

where we have identified S± = (5'j_, 5^, S^) . 

It is the term involving the derivative of delta function in (^^) that leads to difficulty 
in quantization of the system. The proposal of Faddeev and Reshetikhin is, therefore, to 
introduce the new Poisson bracket 

{SUx),Si{y)} = -e''''SUx)S{x-y), 
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(2.5) 

{Sl{x),S'Liy)} = 0, 

which does not contain any non-ultralocal terms and a new Hamiltonian which would give 
rise to the same equations as in (2.4). It follows now from ( p.5| ) that the consistency of 



these relations requires 

= = (2.6) 



On the other hand, for the system to have the same equation of motion ( |2.4| ), the Hamil- 
tonian of the system must modify as well (along with the Poisson bracket) and the new 
Hamiltonian is given by 



n 



- jclxS^- (2.7) 



where Vs^ and Vs_ are the momenta associated with the two variables, with the explicit 
forms: 



J \s[Tsl 



\s\ + s^_ 

It is easy to check using (|2.5D that 

{r,,s+} = -dJ+, {Vg^,S-} = -dJ., (2. 



so that the dynamical equations (2.4) arise as Hamiltonian equations with the modified 
Hamiltonian (2.7) as well as the modified Poisson brackets ( |2.5[ ). 

The dynamical system described by (^^) has an infinite set of conserved charges 
associated with it and the Hamiltonian (|2.7|) can be determined from the low order charges. 



To see how the infinite set of conserved charges arise, let us note that the dynamical 



equations (2.4) can be obtained from the zero curvature (flat) condition associated with a 



one parameter family of currents 

dtS+{X) - dxS^iX) + 25+(A) x S_(A) = 0, (2.10) 
where A is a constant (spectral) parameter and we have defined 

5'±(A) = -^±-^, (2.11) 
A — a A + a 

with a = 1/7. The transfer matrix 

T(x, y,X) =V exp I J dz S^{z, A) 
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where V denotes path ordering, can be decomposed in the standard form 



T{x, y, X) = [I + W{x, y, A)] exp (Z(x, y, A)) [/ + Wix, y, X)]-' , (2.12) 

where Z(a;,y,A)and W{x,y,X) are respectively diagonal and anti-diagonal matrices. The 
monodromy matrix, which contains all the conserved charges of the theory, is defined as 
(we assume that the theory is defined for —L<x<L with the continuum limit obtained 
for L oo) 

Tl{X)=T{L,-L,X) (2.13) 

One can now obtain the local charges by expanding the monodromy matrix around the 
two poles in ( p. 11 ) which leads to two sets of conserved charges 9?^^"^ from the series: 



V'L (A) = ^ (A ± a)>J("\ (2.14) 

n 

where the generating functional (pL (A) is defined through the relation 

Ti [Tl {X)]= 2 cos if l{X). (2.15) 
For the FR model, the decomposition ( 2.14| ) has the explicit form: 



(A) = Y^^t^''^ + + (A ± a) ^t^'^ + ... (2.16) 
In this series the trivial conserved charges correspond to the conservation of total spin: 

L 

±(-1) 



= J \S\dx (2.17) 

-L 

The first non-trivial charges have the form: 

±(0) T. -L ^ f ^+ ' ^ 



where Vs± are the momenta defined in (|27 

To summarize, therefore, the FR proposal is to modify the Poisson bracket (without 
any non-ultralocal term) as well as the Hamiltonian of the theory such that the same 
dynamical equations result. The theory can now be quantized. The idea here is that the 
term with the derivative of the delta function in the original Poisson brackets (say, in ( |2.2| ) ) 
is an anomalous term (can be thought of as an anomaly) which may arise in some limiting 
manner after quantization. 

3. String Sigma Model. 

As we have mentioned earlier, in the Green-Schwarz formulation, the superstring on AdS^ x 
can be described as a symmetric space sigma model [66]. Let g denote an element of 



-6- 



the graded group PSU{2,2\4:) (namely, it represents a map from the worldsheet onto the 
graded group P5'C/(2, 2|4)). Then, the left-invariant current associated with the group 
J = —g~^dg can be decomposed in terms of the superLie algebra elements of psu(2,2|4), 
which has a natural Z4 symmetry (automorphism), as 

J = -g-^dg = H + P + Q^ + Q'^. (3.1) 

Here H denotes elements of the maximal (non-compact) bosonic subalgebra so(4, 1) x so(5), 
while P denotes elements of the bosonic complement and Q^, represent the Grassmann 
elements of the superalgebra under the Z4 grading. In terms of these variables, the string 
sigma model action on AdS^ x can be written as 



^ j str{P A*P -Q^ AQ"^) , (3.2) 
where "str" stands for supertrace, "*" denotes the Hodge dual operation and the second 



«=2 



term in (3^) represents a fermionic Wess-Zumino term. Introducing the one parameter 
family of flat currents [18] j(A) = — ^^^(A)d^(A), where A is a constant spectral parameter. 



i(A) = ^ + l^^+T^*^+ VS ^ + Vr5 ^ ' ^'-'^ 

such that J(A = 0) = J, it is easy to check that the vanishing curvature condition for this 
current 

dJ-JAJ = 0, (3.4) 



leads to the equations of motion for the system following from the action (3.2). Here we 
have defined Q = + Q^,Q' = - Q^. 

In [18,19] we have calculated explicitly the algebra { Ji(fT, t)^Ji{a' , t)}, and have shown 
that this Poisson bracket algebra contains the problematic da6{a — a') term in general. 
Therefore, the FR approach may be relevant in understanding the quantization of this 
system. The presence of the Wess-Zumino term in the action in (^^) can also be incorpo- 
rated into the proposal of FR [32], by weakening the constraint ( ^.6| ) so that the lengths 
1 5+ 1 and are not equal. Thus, the Faddeev-Reshetikhin proposal is a promising and 
interesting scheme that may potentially allow one to proceed with the quantization of the 
string sigma model on AdS^ x and deserves to be investigated. However, to keep the 
discussion parallel to the SU{2) principal chiral model of the last section, we will restrict 
ourselves here to the string sigma model on the i? x subsector of AdS^ x . Unlike the 
principal chiral model in flat space, here the equations of motion need to be supplemented 
by the Virasoro constraints. Essentially, what this means is that, unlike in the case of the 
principal chiral model, in the case of the string, the constraint ( |2.6D does not have to be 
set by hand. 

In fact, an important progress has been made in this direction recently by Klose and 
Zarembo [41], where the S-matrix calculation has been carried out for the FR model (among 
several other models). The calculation involves summing up a particular set of Feynman 
diagrams, the bubble diagrams, for the two-particle scattering, which are the relevant 
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diagrams for calculations performed in the wrong vacuum. Using the factorization property 
of the S-matrix of an integrable system [53], one can then write down the S-matrix for 
general A^-particle scattering. The corresponding Bethe Ansatz for the theory can then 
be obtained by imposing relevant boundary conditions. Since this calculation is rather 
relevant from our point of view, here we briefly review their main results directing the 
reader to [41] for technical details. 

The action for strings propagating on i? x 5"^ subsector of AdS^ x has the form: 



S = -— [ drda qf"" 
Att J 



(3.5) 



where -^/AySTT represents the string tension, the current J = —g ^dg G su{2), and g is an 
element of the standard map S'^ — > SU (2) 



+ iX'^ X^ + iX'^ 
-X^ + iX'^ X^ - iX"^ 



(3.6) 



The action (^]^) is written in the conformal gauge 

^/=5ff'^"=r/'^^ (3.7) 

and introducing the light-cone coordinates cr-t = ^(r it o"), = (?t- ± we note that the 
Virasoro constraints of the theory take the form: 

Tr{Jl) = -2{d±X^f . (3.8) 

Furthermore, if we use the parameterization 

X° = KT, (3.9) 

where k denotes a constant, the Virasoro constraints in ( |3.8D take the simpler form: 

rr(4) = -2k^. (3.10) 

To make connection with the SU{2) principal chiral model considered by Faddeev and 
Reshetikhin, let us identify: 

J_l_ = iKS± ■ a, 

(3.11) 

s± = {s}^, S^, S^), 

so that the Virasoro constraints and the equations of motion of the theory, in terms of the 
spin variables S± become 

Si = 1 

d^S± ^ kS+ X S_ = 0, (3.12) 
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which can be compared with (2^) and ( |2.6| ). Thus, the analogy with the FR model is now 
clear and one can follow the FR proposal to quantize the system. The action corresponding 
to the Hamiltonian (|2.7|) takes the form: 



S 



(fx 



C+{S-) + C^{S+)]--S^-S. 



K -> 



(3.13) 



where C+(S'_) and C-{Sj^) are the Wess-Zumino terms which have the explicit forms 

1 



C.{S, 







£"''^318^318-31, 



(3.14) 



with the boundary conditions 

S±{T,a,C = l) = (3± 



const, 



S±iT,a,C = 0) = 3±iT,a). 



(3.15) 



(3.16) 



Rewriting the Wess-Zumino terms ( 3.14 ) in the local, non-covariant forms leads to the 
expression in (|2.8| ) for the momenta of the theory, if one chooses \3\ = 1. Solving the 
constraint 5^ = 1 and introducing 



3i + i3l 



V2^{1 + 31) 

the action ( 3.13| ) can be written in the unconstrained form: 



(3.17) 



3 



8+<t). - <P-8+^*_) + - {(t^\8-<P+ - + K (|0+|2 + |</>_|2) 



(3.18) 



-K^{l-\<l)+\^){l-\(t>-m<t>\(t>- + (t>-<t^+) - 2k|0+|2|0_|2 
Finally, introducing the two component bosonic spinor 



the action (|3.18D can be recast into the following compact form for terms up to the order 
(f)^ (quartic order): 



3 = Jd^x [i^rD^cl) -m^^-g ((/.7'^<A) {^J^4>) + 0(/)] , 



(3.19) 
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where Dq = do — im — ig(j)4>; D\ = di = dx ] m = k\ g = ^. The non-covariance of the 
model is hidden in the field dependent chemical potential in the definition of the covariant 
derivative, which when written out exphcitly would lead to a term in the quartic interaction 
Hamiltonian that violates Lorentz invariance. We note here that throughout the paper, 
the interactions are assumed to be normal ordered, although we do not explicitly write the 
normal ordering symbol for simplicity. 

We do not go into the details of Klose and Zarembo's calculations which are explained 
nicely in [41]. Rather, we would like to summarize the essential features of their calculation 
which would be relevant for our discussion. First of all, the theory in ( [3.19 ) has both 
positive and negative energy solutions, as any free relativistic theory would have. For free 
particles, they satisfy the momentum space equations 



m)u{k) = = (^ + m)v{k), 



and have the explicit forms 
u{k) -- 



m 



62 

e 2 



v{k) 



m 



— 62 
6 2 



(3.20) 



(3.21) 



Here we have used the rapidity variable defined as 



= msinh/3, = \/ k"^ + m? = m cosh /?, (3.22) 
to parameterize the solutions and we have chosen the Lorentz invariant normalization 

u{k)u{k) =2m = -v{k)v{k), (3.23) 

with u{k) = u^{k)'j^. 

If we quantize the theory in the wrong vacuum: 

(j){x)\0)=0, (3.24) 

then the field can be decomposed completely in terms of annihilation operators as 



4>ix) 



dk 



-ik-x 



u{k)a{k) + e'^-''v{k)h{k)] , 



(3.25) 



Here we have identified k^ = Ej^,k^ = k and using the (nontrivial) equal-time commutation 
relation for the fields 



(3.26) 



it can be checked that the nontrivial commutation relations for the creation and annihilation 
operators take the forms 



a{k),a\k')\ = 5{k - k') = \h{k),h\k') 



(3.27) 



Normally, in a relativistic theory, the time ordered propagator is the Feynman (causal) 
propagator. However, as a result of quantizing the theory in the wrong vacuum, the 
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Figure 1: The bubble diagrams contributing to the S-matrix 



propagator of the (relativistic) theory becomes retarded (This is completely in spirit with 
the Bethe ansatz methods in other systems) and has the form 

D(k) = % . (3.28) 

As a result, the only diagrams that contribute to the two particle scattering matrix are the 
bubble diagrams in Fig. 1 . (This also clarifies why the quartic part of the action is sufficient 
to study scattering of two particles.) This simplification is automatic in non-relativistic 
models where propagators are by definition retarded. However, in relativistic models, this 
is achieved by quantizing the theory in the wrong vacuum (the true vacuum being the one 
where all the negative energy states are filled). 

The calculation of the bubble diagrams is quite simple once the 1-loop diagram is 
known. The interaction vertex in this theory can be written in a tensor notation as 

= ^ [7° ® 1 + 1 ® 7° - 27^" ® 7m] > (3-29) 

and factoring out the vertex parts, the one loop contribution to the bubble diagram has 
the form 

1-loop = j ^ D{p + p'-q)® Dig), (3.30) 

where D{q) is the retarded propagator defined in ( [3.281 ). For positive energy external 
particles, this integral can be evaluated to yield (for technical details, we refer the reader 
to [41]) 

1-loop = „ n . S7V K/ + m) IS) + m) + + m) IS) (/ + m)] , (3.31) 

8m^ smh(p — p'j i-v / \ /j 

where /3 is the rapidity defined to be (see ( 3.22| )): 

= m cosh P; p^ = m sinh /S . (3.32) 

The general n-loop contribution to the S-matrix element , coming from the bubble dia- 
grams, is then obtained by raising the one loop result multiplied by the vertex functions 
raised to the nth power and taking the matrix element between the incoming and the out- 
going positive energy states \p,p') and \k, k'). The complete two particle S-matrix element 
corresponds to summing over all loops and yields the result (because of energy-momentum 
conservation, k = p, k' = p' or k = p',k' = p as explained in [41]) 

{k,k'\S\p,p') = Si(5,(5') = ^±^, (3.33) 
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where 

cosh(^) 

^ ~ " '/3-/3' 



cosh( ^'^^ ) — cosh( ^ ^ 
sinh(^— 2- 



A^g \ [,Z7? ' • (3-34) 



Assuming the integrabihty of the model, and using the factorization property of the S- 
matrix [53], one can now write a general N -particle scattering S-matrix as a product of 
two-particle S-matrices. It is worth noting here that the calculation of the S-matrix el- 
ement involves an inner product between the initial and the final positive energy states. 
As a result, if the interactions of the incoming positive energy particles generate an inter- 
mediate state that is orthogonal to the positive energy states, the inner product with the 
positive energy outgoing states will not see it. Therefore, the S-matrix calculation will not 
be sensitive to such an issue. Normally, this is not an issue in other integrable models. 
However, as we will see in the next section, this question becomes quite important in the 
diagonalization of the present Hamiltonian. 

Once the two particle S-matrix is known, one can write the Bethe equations as: 

As we discussed in the introduction, the next logical step to carry out is to identify the 
negative energy states and fill the Dirac sea with the purpose of constructing the physical 
vacuum and the physical S-matrix. This well-known procedure has been worked out in 
detail by Korepin for the fermionic Thirring model [57,67]. As the first step towards this, 
one needs to know the explicit form of the two particle wave- function which will diagonalize 
the Hamiltonian. This is the question that we take up in the next section. 



4. Diagonalization of the Hamiltonian 

Before proceeding, let us fix our conventions and notations. We have already mentioned 
that we use the Bjorken-Drell metric, ry^'^ = (+, — ) and although the explicit forms of 
the Dirac matrices are not relevant, for completeness we note that we use the following 
representations for the 7 matrices 

l^ = (yi = i^i^^ ' 7^ = -ic^2= ,75=7V = c^3= ■ (4.1) 

The two component spinor indices, when used, will be labelled by the beginning of the 
Greek alphabet a, /?, • • • = 1, 2. We use natural units where h = c = 1. 

The Hamiltonian density following from the quartic action in (|3.19| ) has the form ^ 

TL = TLq + TLi 

= -i<l)^a3d,^ + m<l)^ai^ - g (</.tai</>) - ((AV)(</'V) + ((AV3</.)(</.V30)1 , (4.2) 



^As is conventionally done, we have omitted the term m(cf>^(f>), which is proportional to the number 
operator, and commutes with the Hamiltonian. 
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which leads to the Hamiltonian 

H = Ho + Hi= I dxH. (4.3) 



As we have noted earher, this is the Hamiltonian that is responsible for the dynamics in 
the two particle sector. We note, that the structure of the Hamiltonian ([4. 21) is reminiscent 



of the bosonic Thirring model, namely, the last two terms in (4^) describe exactly the 
interaction of the bosonic Thirring model (while the first term in the interaction Hamil- 
tonian violates Lorentz invariance). Thus, we can follow the standard procedure for the 
diagonalization of this Hamiltonian [47,51]. 

The theory, like all relativistic theories, has positive and negative energy solutions and 
let us construct the one particle states of the theory. As described in the last section, we 
will assume that the vacuum is annihilated by the field operator, 

(j){x)\{)) =Q. (4.4) 

Let us write the positive energy single particle state with momentum k as 

= j dx xL+)(x|A:)4(x)|0). (4.5) 

Requiring this state to be an eigenfunction of the Hamiltonian ( [4. 3D with positive energy, 

H\m)[+)=Ek\m){+), (4.6) 

determines the spinor x{x\^) to be (basically here one commutes the field variables to the 
right using ( ^.26 ) until it annihilates the vacuum ( 3.24| )) 



/3 



X^+\x\k) = x^+^(x|/?) = n+(A:)e*'=^ = % gi(msinh/3)x^ ^^^^^ 



e 2 



where we have used the rapidity variable defined earlier in ( |3.22 ), namely. 



= msinh/3, Ej^ = \/ k'^ + m? = m cosh (3, (4-8) 

to parameterize the solution. On the other hand, for the negative energy single particle 
state with momentum k, defining 

= j dx xViAmiixm, (4.9) 

and requiring that 

iJ|V'(A;))(^) = -£;fe|V(A:))(_), (4.10) 

we determine 



X^-\x\k) = x^~\x\(i) = u.{k)e'^'' = ^/^i V ] e'^~'^'^^\ (4.11) 
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The spinors in (|4.7| ) and ( 4.11 ) are normalized to be consistent with the invariant 
normahzation in ( |3.23| ) and, in fact, we can relate them to the covariant spinors in dPI) and 
(3.21) simply as u+{k) = u{k),u-{k) = v{—k). Thus, these are, in fact, free single particle 
states, although they are eigenstates of the full Hamiltonian. This can be understood from 
the fact that the interaction Hamiltonian leads to a vanishing contribution to the single 
particle eigenvalue equation. We note now that if we identify the annihilation operators 
associated with the original field variables as 

dx 



aa{k) = 

then using ( 3.26| ), we can obtain 



e-'^^M^), a = 1,2, 



-ikx 



a,P= 1,2. 



(4.12) 



(4.13) 



Furthermore, let us define a new set of annihilation operators through a change of basis as 



A2{k) 



1 



1 



^{k)aa{k) = cos6kai{k) +sm9ka2{k) 



ul_ ^{k)aa{k) = — sin 0kai{k) + cos 6*^02 (A;) 



dx 



-ikx A 



./me; 

dx 



V^^k 



ul^{k)4>{x), 

(4.14) 

*^^nl(A:)0(x), 



with 0fc satisfying the relation ktaii26f; = m. In fact, from the definition of the positive 
and the negative energy spinors in (4/7) and ( 4.11 ), we note that 



62 



e 2 



cos Ok 



V2cosh/^' 



sm 



(4.15) 



V2 cosh 

In this new basis, it is easy to check that the free Hamiltonian is diagonal with the form 

dk Ek {A\{k)Ai{k) - A\{k)A2{k)) , (4.16) 



(4.17) 



Ho-- 

and we can identify the single particle states in ( [4.5[ ) and ( [4.9|) with 

4(A:)|0) = -^|V;(A;))(+), 4(fe)|0) = -^|V;(A;))(_) 



Namely, ^|(A;) and ^2(^) '"^'^ thought of as the creation operators for single particle 
states with positive and negative energy respectively. It is worth noting here that the 
relations ( [4.14] ) can also be inverted to give 



ai{k) = cos0kAi{k) — sm9kA2{k), 
a2{k) = sm9kAi{k) + cos9kA2{k). 



(4.18) 



More importantly, using the structures of the positive and the negative energy spinors, we 
can also invert (4.14) to determine the field expansion in terms of the new operators as 



<A(x) 



dk 



Akx 



{u+{k)Ai{k)+u-{k)A2{k)). 



(4.19) 
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To proceed with the diagonahzation of the Hamiltonian in the two particle sector, we 
follow the standard procedure. However, to see clearly the effect of the interaction term in 
( |4.2| ) which violates Lorentz invariance, we put an arbitrary parameter in front of it and 
write 

Hi = -g / dx (a0V0^o-i<A - '/'V'/'V + 0-3 Vac/)) , (4.20) 



where we note that for a = 1 we have the quartic Hamiltonian of the FR model, while 
for a = we get the bosonic Thirring model. Using the relations in ( 4.14| ) as well as the 
definition of the single particle states in ( [4.17] ), we note that in the absence of interactions, 
the two particle positive energy state with momenta fci, k2 can be written as 

\ki,k2)(^+) = 4(A:i)4(fc2)|0) 

(4.21) 

In the presence of the quartic interactions, however, these would no longer correspond 
to the two particle (positive energy) eigenstates of the complete Hamiltonian. Thus, we 
generalize the definition of the two particle state with the ansatz 

miM))(+) = I -^^^^x(xi,:r2|A:i,fc2)e*(^'^^^+'^^^^)(</.t(xi)n+(A;i))(</.t(x2)tx+(A;2))|0), 

J ATTy/Ek-^Ek2 

(4.22) 

where we assume 

x{xi,X2\ki, k2) = 1 + iXe{xi - X2). (4.23) 

Here A is a function of the momenta as well as the interaction strength (such that it 
vanishes in the absence of interactions) to be determined. One can, in principle, take a more 
general modification of the wavefunction, which we have done. However, the difficulty in 
diagonalizing persists nevertheless and, therefore, we discuss the issue with the conventional 
form of the generalization for the state. 

The calculation, which is slightly tedious, can be carried out in two steps. First, using 
(3.26) and moving the field operators to the right until they annihilate the vacuum ( 3.24| ), 
we obtain 

Holtpiki, A:2))(+) = {Ek, + Ek,)\il>{ki,k2))(+) + |^), (4.24) 

where 

\R) = _ 4Asm( gfc,-gfc , ) I ei(fci+A.).<^t(^)<^t(^)|o). (4.25) 



An^Ek.Ek^ J 

In carrying out the calculation, we have used the standard regularization 5{x)e{x) = 0. 
The action of the interaction Hamiltonian in ( [4.20 ) can similarly be calculated and leads 
to 

Hi\^P{k^,k2))M = -^^!i^^i|L±£^ f dx e^^^-^^-> (0t(x)0t(x) + 4(x)4(x)) IC^.26) 



AiT^Ek^Ek^ 
2g (g cos(6'fc^ -Ok^) -2 sin(6'fci + 6^^ J ) 
An^Ek^Ek^ 
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The two particle state ( [4.22| ) would clearly be an eigenstate of the full Hamiltonian only if 
the sum of the terms in ( |4.25D and (4.26) cancel, which would also determine the parameter 
A as a function of the momenta and the interaction strength. However, in the present case, 
we see that since the structures of the terms in the two expressions are quite different, such 
a cancellation is not possible unless a = which would correspond to the bosonic Thirring 
model. In this case, we can determine 



A = <7 



sm 



-g coth 



(4.27) 



where, using (4.15), we have expressed the result in the rapidity variables in the last form. 
The S-matrix can now be written as 



S 



1-iA' 



(4.28) 



which is well known for the bosonic Thirring model. This result also holds for the fermionic 
Thirring model [68,69]. However, we also note that for any nontrivial value of the Lorentz 
violating parameter a, the two expressions in ( [4.25| ) and ( 4.26| ) cannot be cancelled and 
hence the Hamiltonian cannot be diagonalized in the two particle sector. 

To make connection with the field theoretic calculations in [41], let us observe the 
following. We can rewrite the terms on the right hand side of ( [4. 261) as 



dx 



J{ki+k2)x 



47r/B^ 



k2 



([extra) - ig{acos{ek, - ^fcj " M^k^ + ^fej).^! (x)4(x)[0)) , (4.29) 



where 



[extra) = -ga (sin(6lfcj + 6lfc2)(4(x)0| (x) + 4(2;)4(^) - 2cos(6'fc, - 6lfc2)0| (x)4(2;) j [0). 

(4.30) 

Using the definitions in (|]2l|), ( pT^ ) and ( pT^ , it is now straight forward to check that 



(+){ki, A;2 [extra) = 0. 



(4.31) 



Therefore, if we take the inner product of the sum of ( |42^ ) and ( ^^261 ) with a positive 
energy two particle state, the sum will vanish provided (we remind the reader again that 
because of energy-momentum conservation, the momenta of the out states will coincide 
with a permutation of the momenta of the in states, as is well known [41]) 



A 



a cos (6*^1 - 9k2) - sm{9k^ + 6*^2 



acosh^i|^ 



cosh^i^ 



sm 



^k2) 



sinh 



/3i-/32 



(4.32) 



which, for a = 1 (FR model) reduces to the field theoretic result in ( 3.34| ). This is indeed 
what we have tried to point out in the introduction. Namely, the S-matrix calculation 
involves calculating matrix elements and, consequently, is not sensitive to states orthogonal 
to the external states generated in the intermediate steps. These states, on the other hand, 
may be very important in the study of the diagonalization of the Hamiltonian. We will 
investigate this question in more detail in the next section. 
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5. Solving the Boundary Condition 



Although the analysis of the previous section makes it clear that the Hamiltonian is not 
diagonahzable in the presence of Lorentz violating terms, neither the origin of the problem 
nor the possible remedy is very clear. For that reason, let us analyze the reason for the 
difficulty in diagonalization from a different point of view. In this section, we will investigate 
the quantum mechanical problem associated with the diagonalization of the Hamiltonian 
in the two particle sector. 

Let us write the quartic Hamiltonian for the system ( [4 .2] ) (with the interaction in 
( |4.20D ) in the form 



H 



j dx (^-i(f>^ a3dx(l) + m(j)^a^(j) - gVa/3,^s (t>i(pl4>-y4>s^ , (5.1) 



where V is the tensor related to G defined in ( ^.29 ) and has the explicit form (Unfortunately, 
the Lorentz violating parameter is called a just like the spinor index simply because we 
are running out of letters. However, we believe that there will be no confusion because of 
this.): 



/O 

a 
a 



a 
-2 
-2 

a 



a 0\ 
-2 a 
-2 a 

a J 



(5.2) 



We will use a tensor product notation in this section which makes the difficulty associated 
with the problem of diagonalization more transparent. Thus, the outer product of two 
positive energy spinors ( |4.71 ) will be represented in this notation as a 4- component column 
vector 

i 2 \ 

/3l-/32 

C 

U++,al3{kl,k2) = U+^a{kl)u+^j3{k2) = m 



01-02 
2 

01+02 
2 



(5.3) 



(and so on for other spinors) and the contraction of the spinor indices would simply corre- 
spond to matrix products in this notation, which simplifies the calculations enormously. 
Let us take a general ansatz for the two particle positive energy state as 



|V'(A;i, fc2))(+) = / dxidx2 Xai3{xi,X2\ki,k2)(t>i{xi)(t>p{x2)\0), 



(5.4) 



where, unlike in ( 4.22| ), we have left the form of the wavefunction Xa/si^i, X2\ki, k2) arbi- 
trary, to be determined from the equations. Requiring that this state represents the two 
particle state of the complete Hamiltonian with positive energy of the form 

H\i,{ki,k2))^+) = {Ek,+Eu,M{kiM))i,+). (5.5) 

leads to (once again, one simply commutes the field variables to the right using ( 3.26| ) until 
they annihilate the vacuum ( 3.24| )) 



dxidx2 



i-i ((0-3 l)«9xi + (1 ® cr3)dx2) +m{ai^l + l^ai) 



-17- 



-{Ek^ +Ek2){l (g) 1) - 2gV6{xi - X2))^f^^^s Xjsixi, X2\ki, (j)i{xi)<l)^p{x2)\0) 

= 0. (5.6) 

Thus, we can interpret the expression in the large parenthesis (without the energy eigen- 
value terms) as the quantum mechanical Hamiltonian Hqm in the two particle sector. To 
determine the two particle state (^.4[), we need to solve the quantum mechanical equation 



-i (((73 O l)5^j + (1 O 0-3)5^2) + m (fJi 1 + 1 O cTi) 



-iEk,+Ek,){l 



— X2)]a(3,'y5 X'y5{xi-,X2\ki,k2) 



0. 



(5.7) 



Equation ( ^.71) is a first order equation with a delta potential whose coefficient (strength) 
has a rather nontrivial tensor structure. In this case, we expect that the wave function itself 
will be discontinuous across the boundary xi = X2- Away from the boundary (namely, in 
the region xi < X2 or xi > X2), the wave functions will correspond to the free two particle 
positive energy solutions whose coefficients must be determined by matching the disconti- 
nuity across the boundary. Thus, let us choose the conventional general wave function of 
the form^ [51]: 



Xap{xuX2\ki,k2) = e*(^i"i+'=2"2)(l + iAe(xi-X2))C/++,a/3(A:i,A:2) 



(5.8) 



where C/++^q/3(A;i, ^2) is the outer product of two positive energy solutions defined in (^). 
Here A is assumed to be space independent and should, in principle, be determined from 
matching the boundary condition. We comment here that the wave function in (|5.8| ) can 
be made even more general by giving A a nontrivial tensor structure. However, as we would 
see shortly, this does not help in the matching of the discontinuity across the boundary. 
It follows from (|5.7|) that at xi = X2, the discontinuity has to satisfy 



(5.9) 



((73 (g) 1 - 1 Cr3)al3,'y5 (2^) {U^+^jsiki, ^2) - U++^ys{k2, ki)) 
= '^9Va(s,^5 (^7++,75(A:i, k2) + U++^^s{k2, ki)) . 



As a result, the discontinuity relation takes the explicit form (in the outer product notation 
introduced in (K.3D), 



4A 



/ 
sinh 
sinh 



/3i-/32 



V 







25 



a cosh ^^"2^ 



\ 



V 



a cosh Ml _ 
a cosh ^i±^ - 
a cosh 



2 cosh 



- 2 cosh 



(5.10) 



Since the left and the right hand sides of the equation have quite different matrix structures, 
it is clear that the discontinuity relation cannot be satisfied and, therefore, a solution to 



For simplicity, we omit here the multiphcative factor of — 
analysis. 



which is not relevant for our 
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( |5.7| ) cannot be obtained for a 7^ 0. For a = 0, for which the model corresponds to the 
bosonic Thirring model, the discontinuity condition can be satisfied and determines 

/5i-/?2 



A 



-g coth ■ 



(5.11) 



in agreement with the result (4.27) in the operator method. 

At this point, one may wonder as to whether one cannot take a more general ansatz 
for the wave function in ( |5.8| ), for example, by giving a tensor structure to A allowing for 
a more compatible structure on the left hand side of (5.10). To investigate this question, 
let us note here that 



(T3 ' 



0-3 



/O 

2 



\0 



o\ 


-2 
0/ 



(5.12) 



is a very special operator which, acting on a four component column vector, removes the 
top and the bottom elements. Therefore, even if we give a nontrivial tensor structure 
to A to allow for a more general component structure for the spinors, when the operator 
((T3 1 — 100-3) acts on it, it would project out the top and the bottom elements and bring 
the spinor to the form on the left hand side of (|5.10| ). Therefore, matching the discontinuity 
relation is the main problem because of which the Hamiltonian is not diagonalizable. 

It is worth reflecting here on the connection between our analysis and the calculation of 
the S-matrix [41]. Let us note that although the discontinuity relation for the wave function 
(|5.9|) cannot be matched, if we take the inner product with the two particle positive energy 



state U 



ii,k2), equation ( [5.9D can determine the parameter A to be 

cosh^i^ 



A 



9 



acosh^i|^ 



sinh 



/3i-/32 



(5.13) 



which is what we have obtained in ( 4.32 ) and which, for a = 1, reduces exactly to the 
value calculated in ( 3.34| ) from field theoretic methods. To understand this better, let us 
note that the projection operator for the two particle positive energy states with momenta 
ki, k2 (as in (|5.3|) ) can be easily computed to have the form 

e-/32 
^ -ft 



P++{ki,k2) 



e e 

-{131+132) g 



1 

(ft -ft) 

f3i 



e'- 

e(ft-ft) 
1 

p-ft 



3(/3i+/32) \ 



jft 

af32 

1 



(5.14) 



so that 



-{ki,k2)U++{ki,k2) = U++{ki,h 

With this, we note that the column vector on the right hand side of (| 
decomposed into the sum 



45 







acosh^i±^ 
a cosh Mi 



cosh 



ft -ft 



cosh^i^ 



/ cosh^i^ \ 



+ 2ag 



cosh 



/3i+ft 



■ cosh Mi 

ft -ft 



y cosh 



(5.15) 
can be uniquely 



(5.16) 
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where the second vector is annihilated by the projection operator i-*f+ in ( ^.14| ). In the 
absence of this second term, the discontinuity relation ( ^.10 ) (or ( |5.9|) ) can be solved and 
yield ( 5.13| ). This is the reason why taking the inner product with positive energy states 
allows us to solve the discontinuity relation leading to the result from the field theoretic 
calculation. The important thing to note here is that in the perturbative calculation, one 
is evaluating matrix elements between positive energy states and, therefore, the calculation 
will not be sensitive to state vectors that are orthogonal to such states if they are generated 
in the intermediate steps. 

There is a second way to look at this issue which is quite interesting. Let us note that 
we can decompose the interaction potential (|5.2|) uniquely as 



V = y + 



(5.17) 



where 



V 



y(o) 



a 



/o 












o\ 


a 




-1 




1 


a 


a 




-1 




1 


a 


Vo 












y 


/ 





1 


1 








1 










-1 




1 










-1 


V 





1 


1 




0/ 



(5.18) 



These matrices have the property that 



VU-^.-^-{kl, = 4m(acosh 



P1 + P2 



cosh 



/3i-/32 



1 
1 



y(°)[/++(A:i,A:2) = 2ma 



I 



\ 



cosh^\ 
cosh^ 
cosh Mi 
cosh^i^ 



(5.19) 



This is precisely the decomposition of the state vectors that we have discussed in ( 5.16| ). 
However, here the decomposition is in terms of the potential. What we see is that 
acting on a two particle positive energy state would lead to a state that is orthogonal to 
such a state. As a result, the S-matrix calculation, whether it is carried out with the vertex 
involving V or the full vertex involving V ^ would lead to the same result since the extra 
terms generated by V^^'^ would drop out in the matrix element between positive energy 
states. On the other hand, from the point of view of diagonalizability, it is only V that 
can be diagonalized and not the full V . We note that the extra term that is orthogonal 
to the positive energy states is proportional to the Lorentz violating parameter. Such a 
term is not present in the conventional (relativistic) integrable systems and this is a new 
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feature in this model. This is the basis of our claim in the introduction that while a 
diagonalizable Hamiltonian leads to the S-matrix, having the S-matrix (say, from a field 
theoretic calculation) does not automatically imply that the Hamiltonian is diagonalizable. 



6. Diagonalizability and PT Symmetry 

The analysis of the last section is quite interesting and leads to the natural question as to 
what is the most general quartic Hamiltonian within this context that can be diagonalized 
and what is the corresponding S-matrix. Such an anlaysis would also determine the inter- 
action Hamiltonian (potential) that would lead to the S-matrix calculated by Klose and 
Zarembo. This can be carried out systematically along the lines of discussion in the last 
section and we find that the most general Hamiltonian that can be diagonalized has the 
form 

H = Hq + Hi = I dx 



Here a, (3 and 7 are arbitrary real parameters (for the S-matrix to be unitary). We note 
that if a = 7 = and (3 = —1, this model reduces to the bosonic Thirring model that 
has been studied extensively. For a = P = and 7 = — 1, this model corresponds to 
the bosonic chiral Gross-Neveu model which, to the best of our knowledge, has not been 
studied in the literature. Both these models involve Lorentz invariant interactions. Finally, 
the term with the parameter a clearly breaks Lorentz invariance, but does not coincide 
with the Lorentz violating term in (f4.2[l| ), rather generalizes it. 

The quantum mechanical potential, in this case, can be worked out in a straight forward 
manner (see last section) and has the form 

/O o\ 

a 2/? 27 a 
a 2/? 27 a 
\0 0/ 

Substituting this potential into ( 5.10| ), it is clear that the discontinuity relation can be 



V 



(6.2) 



satisfied with 

a cosh ^ „ ' y,^ • 1 1 — " .) 
A = 5 ^ . —■ (6.3) 



a cosh Ml + (/3 + 7) cosh 
sinh^i^ 



2 

It is important to note from ( |6.3| ) that A is an antisymmetric function under the exchange 
of momenta, which is consistent for the Bose symmetry of the wave function. The S-matrix 
for this general diagonalizable model, then, follows from the standard relation 

5=1±^. (6.4) 

It is now clear that the parameters a, /3, 7 must be real so that A is real and correspondingly 
the S-matrix is unitary. 

There are various special cases that one can study from (|6.3| ) and we list only three 
that we think are interesting. 
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1. First, if a = and /? + 7 = — 1, 



\ = -g coth ^ , (6.5) 



which we recognize as the A for the bosonic Thirring model ( [4.27 ). However, what 



is striking here is that there is a one parameter family of interactions with a = 
and /3 + 7 = —1 which share the same value of A and, therefore, the S-matrix. In 
particular, we note that for [5 = —1,7 = 0, the model corresponds to the bosonic 
Thirring model. On the other hand, for /3 = 0, 7 = —1, the model corresponds to the 
bosonic chiral Gross-Neveu model and we find that both the bosonic Thirring model 
as well as the bosonic chiral Gross-Neveu model (along with the one parameter family 
of interactions) share the same S-matrix element. To the best of our knowledge, this 
has not been recognized earlier. We also note that for a = 0, if /? + 7 denotes an 
arbitrary constant (not equal to unity), this simply scales the coupling constant. 
However, for a = 0, if we also have (3 = —7 7^ 0, the S-matrix is trivial in spite of 
the fact that it is apparently an interacting theory. All of this can be understood 
as follows. In the outer product space of 2 x 2 matrices, the completeness relation 
(Fierz identity) takes the form 

^aP^'^S = 2 [^a&5^ji + (o-a)a<5('7a)7/3] , a = 1, 2, 3. (6.6) 

Using our convention for the gamma matrices in (^]^), we can also write this as 

^afS^^i = ^ [^a&^^p + {l^)a5{lti)^l3 + (Ts)^^ (75)7/3] , /i = 0, 1. (6.7) 



Contracting ( |6.7| ) with the bosonic fields (pa'P'y't'fS't'S under the normal ordering sign, 
we obtain 

4>(j)4>(j) - clyy5(j)4)-f5(p = (i)-i^'(t)(t)"f^(t), (6.8) 

which shows that in 1 + 1 dimension, the bosonic Thirring interaction is equivalent to 
the bosonic chiral Gross-Neveu interaction, which explains the results noted above^. 
As a result of this equivalence, the last two terms in ( |6.lD can be combined into one. 
However, we keep them separate in the following discussion just for completeness. 

2. We note next that for a = 1 and /3 -|- 7 = —1, 

cosh _ cosh , , 



This is exactly the A in ( 3.34 ) (and, therefore, the S-matrix) that has been calcu- 
lated by Klose and Zarembo. We find that, although the FR Hamiltonian is not 
diagonalizable in the two particle sector, there exists a generalized interaction violat- 
ing Lorentz invariance that can be diagonalized and leads to the same perturbative 



^For the massive fermionic Thirring model, it is well known that : -ip^^ ip-ip'y : ~ : Tpipipip 
i/'75'0i/'75'i/' :, which arises from the nilpotency of the fermionic fields. 
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S-matrix as in [41]. In fact, as we have noted earlier, one can add any multiple of 
V^'^^ defined in ( |5.18| ) to this potential which would lead perturbatively to the same 
S-matrix element, but such Hamiltonians cannot be diagonalized. We note here that 
for a arbitrary with /? + 7 = — 1, ( |6.3D coincides with ( [4.32| ). 



3. Finally, we note that if q = 1 and /3 + 7 = 0, 

cosh^i±^ 
sinh 



^ = 9 . , (6.10) 



The integrability of this model, to the best of our knowledge, has not been studied 
earlier. 

To summarize the results of the analysis of this section, we have determined the most 
general quartic Hamiltonian involving scalar fields that can be diagonalized. This involves 
three real parameters a,/3,7 (actually two parameters if we use the equivalence in ( |6.8| )). 
The spectrum of N particle states in such a system is real and is given by 

N 

E = m'^cosh.Pi, (6.11) 

i=l 

where the rapidities satisfy the Bethe equation: 

g^sinh;3,L^-Q^(^.^^.)^ (6.12) 



and the S-matrix (|6.4D , determined from ( |6.3D , is unitary. However, if we look at the 
Lagrangian density or the Hamiltonian of the system, we find that it is not Hermitian. 
For example, we note that the Lagrangian density for the system has the form (we can, 
in principle, absorb the chiral Gross-Neveu interaction into the Thirring interaction using 
(|6.8|), but we keep them separate for completeness) 



£ = -m)4> + g [a07°7^,^07/,</. + (34>j^'(t>hf,(t> + -f{^(t>4>(l) - h54>(f>l5(p)] , (6.13) 

while the Hermitian conjugate is given by (up to a total derivative) 
£t = ^i^f'd^ - m)(l) + a^J°{J^)Uh^.^ + Ph^(l}hi.<t) + - 075^750) , (6.14) 

where we have used the reality of the parameters a, (5, 7. Thus, we see that the Lagrangian 
density is not Hermitian because of the interaction term violating Lorentz invariance. 

It is, therefore, surprising that the spectrum of the theory is real and the S-matrix is 
unitary. This can be understood from the fact that even though the theory (Lagrangian or 
the Hamiltonian) is not Hermitian, it is PT symmetric. Such theories have been studied 
quite a lot in recent years from a variety of points of view [60-64]. Here we describe its 
relevance within the context of this integrable model. First, we note that the simple two 
dimensional quantum mechanical example that is discussed extensively within the context 
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of PT symmetry is reminiscent of our positive and negative energy single particle solutions 
in (4.7) and (4.11). Under parity transformation, P, we note that 

-X, t^t, -k, E ^ E. (6.15) 

Parity is a linear operation and we can define its action on the field space by the relation 

P(t){x, t)P~^ = r]p'y^(j){-x, t), P4>{x, t)p-^ = r?p0(-x, i)7°. (6.16) 

Here r/p is a phase denoting the intrinsic parity of the field. Time reversal, T, on the other 
hand is an antilinear operation defined on the coordinates by 

x^x, -t, k^ -k, E ^ E. (6.17) 

In the field space, the transformation can be described through the action 

TcPix, t)T-^ = r?TC750(x, -t), T^{x, t)T-^ = r]*T^{x, -t)-i^C-\ (6.18) 

where rjT is a phase and C denotes the (Dirac charge conjugation) matrix satisfying 

C-S'^C = -(t'^)^. (6.19) 

Although it is not necessary, we can choose the representation 

C = -^7^ so that C757° = 1. (6.20) 

With these transformations, it is easy to verify that under the parity transformation, the 
Lagrangian density is invariant, namely, 

C{x,t) ^ C{-x,t). (6.21) 

On the other hand, since T is an antilinear transformation, the Lagrangian density is T 
invariant only for real parameters a,fi,^ (which is the case we are considering), namely, 

£{x,t) ^ £{x,-t), (6.22) 

only if a = a* , P = /3*,7 = 7*. Thus, for real parameters, we see that the theory is PT 
symmetric although it is not Hermitian. 

Let us next look at the behavior of the wave functions under this symmetry. We 
note that under the combined PT transformation (see ( |6.15| ) , ( |6.17| ) and remember that T 
denotes an antilinear transformation) 

^ikx ^ g*fex_ 23) 

Similarly, if we choose the phase factors to be unity, namely t]p = rjT = 1; under this 
combined operation 

u±{k) C757%±(A;) = lu±{k) = u±{k), (6.24) 

where we have used ( 0.201) . As a result, we see that the single particle positive and negative 
energy wave functions for the system in ( [4. 7] ) and (4.11) are PT symmetric. Similarly, using 
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the relations ( |6.23| ) and ( |6.24| ) in (^]^) , we find that the complete two particle wave function 
is also PT symmetric (remember the anti-symmetry of A in the momenta). In other words, 
not only is the Hamiltonian of the theory PT symmetric, but so are the wave functions 
of the theory. This implies that the theory is in an unbroken PT symmetry phase which 
guarantees the reality of the spectrum as well as the unitarity of the S-matrix. This is 
indeed a novel demonstration of the relevance of PT symmetry in an integrable system. 

7. Conclusion 

In this paper, we have studied in detail the Faddeev-Reshetikhin model, which is relevant in 
the quantization of strings on AdS^ x S^, in the two particle sector. Although the S-matrix 
of the theory has been calculated using field theoretic methods, diagonalization of the 
Hamiltonian is essential to carry out the Bethe ansatz analysis. We find that the quartic 
Hamiltonian for this model is not diagonalizable in the two particle sector. We show this 
both in the operator description as well as in the description of the underlying quantum 
mechanical system. We trace the difficulty to the fact that the term in the interaction 
violating Lorentz invariance requires a discontinuity in the wave function that cannot be 
satisfied. On the other hand, if one takes the inner product of the discontinuity relation 
with positive energy states, the problematic term disappears leading to the correct S-matrix 
element calculated earlier. Further investigation shows that the interaction (potential) 
generates intermediate states that are orthogonal to the positive energy out states and, 
therefore, cannot be observed in the S-matrix calculation (which involves calculating matrix 
elements), but are quite relevant in the diagonalization of the system. To the best of our 
knowledge, this is a new feature that has not been observed earlier in the study of integrable 
systems. It follows, therefore, that while the diagonalization of a Hamiltonian leads to the 
S-matrix of the theory, the knowledge of the S-matrix element (from, say, a field theoretic 
calculation) does not automatically guarantee the diagonalizability of the Hamiltonian of 
the system. We determine the most general Hamiltonian with quartic interactions that 
can be diagonalized as well as the associated S-matrix. Among various special cases, it 
also includes a generalized Hamiltonian that can be diagonalized with the S-matrix as 
calculated by Klose and Zarembo. We show that although this general Hamiltonian leads 
to a real spectrum and a unitary S-matrix, it is not Hermitian. However, we demonstrate 
that the theory is PT symmetric and that wave functions are also invariant under PT. As 
a result, the theory is in an unbroken phase of PT symmetry which is the reason for the 
reality of the spectrum as well as the unitarity of the S-matrix. 
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